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DIS and the effects of fluctuations: a momentum space analysis 
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Among the dipole models of deep inelastic scattering at small values of the Bjorken variable 
X, one has been recently proposed which relates the virtual photon-proton cross section to the 
dipole-proton forward scattering amplitude in momentum space. The latter is parametrized by an 
expression which interpolates between its behavior at saturation and the travelling wave, ultravio- 
let, amplitudes predicted by perturbative QCD from the Balitsky-Kovchegov equation. Inspired by 
recent developments in coordinate space, we use this model to parametrize the proton structure 
function and confront it to HERA data on ep deep inelastic scattering. Both event-by-event and 
the physical amplitudes are considered, the latter used to investigate the effect of gluon number 
fluctuations, beyond the mean-field approximation. We conclude that fluctuations are not present 
in DIS at HERA energies. 



I. INTRODUCTION 



The well known correspondence between the evolution in QCD at small-x (at high energy) and a reaction-diffusion 
process has been the main source of the recent knowledge concerning the scattering amplitudes and their evolution 
towards the high energy limit. In particular, it has been realized that the Balitsky-JIMWLK equations do not take 
into account the influence of the gluon number fluctuations [li l2) S Ij]. In the mean field approximation, these 
Qh| equations reduce to a unique equation, the so-called Balitsky-Kovchegov (BK) equation ^, [Bj , the simplest equation 
^ • that describes the evolution of the dipole scattering amplitude with rapidity Y = Inl/x. Being a mean field version 
of Balitsky-JIMWLK equations, BK equation does not include the effects of discreteness and consequently of the 
fiuctuations. Among the consequences of the fluctuations, at least in the fixed coupling case, one can cite, for 
example, the slowing down of the approach towards the unitarity limit, as compared to the mean field framework, 
^ and the break down of the geometric scaling [7], a phenomcnological feature observed at the DESY ep collider HERA, 

in the measurements of inclusive "/*p scattering, which is naturally explained in terms of the traveling wave solutions 
of BK equation [S] . 

At small- a;, electron-proton deep inelastic scattering (DIS) can be seen in a particular frame, called dipole frame, 
which allows the factorization of the virtual photon-proton cross section. In this frame, the proton carries most of the 
total energy, but the photon has enough energy to split into a quark-antiquark pair, or a dipole. This qq pair then 
QQ ' interacts with the proton. The virtual photon-proton cross section can be written as 

O- 



•r^. all{Q^Y) = J d'rj dz\^TAr,z;Q^)\ <Jiip{r,Y), (1) 

Cq ' where the labels T and L refer, respectively, to the transverse and longitudinal parts of the cross section, r ^ x — y 
is the vector which gives the transverse size of the dipole, x and y being the transverse coordinates of the quark 
and the antiquark; z is the fraction of the momentum of the photon carried by the quark and '^T,Li'r', z; Q^) are the 
transverse and longitudinal wave functions for the photon to go into the dipole, whose explicit forms are well known 
from QED. Using Eq.([T]) one can obtain the expression for the F2 proton structure function through the formula 

F^ix, Q2) ^ Ql_ L7-P(^, g2) ^ ^7*P(^^ q2)1 (2) 

The quantity Udip in Eq.lU]) is the dipole-proton cross section which can be expressed as 
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{r,Y) = 2 cfb{T{r,b))Y. (3) 



""dtp V ' -^ / — " / "- " \^ ^' , ";/y , 

where {T{r,b))Y is the average scattering amphtude - the notation (•) denotes the average over aU the reahzations 
of the target (proton) color field ~ for the dipole-proton scattering at a given impact parameter b — {x + y)/2. 

As the explicit forms of the virtual photon wave functions are well known, one is left with the parametrization of 
the scattering amplitude, that is, one has to model the dipole-proton cross section. Different approaches have already 
been proven successful, giving good description of the data, among them one can cite the model proposed by Golec- 
Biernat and Wusthoff, the GBW saturation model [£], and the model by lancu, Itakura and Munier, the IIM model 
or CGC fit [lO|, both developed in coordinate space. These models have been recently considered in order to study 
the effects of fiuctuations on DIS in the fixed coupling case [11[ . Although the description of DIS data is improved 
once gluon number fluctuations are included, it is not possible to state that the improvements come from them, but 
they may come from the geometric scaling violation present, for example, in the diffusion part of BK solution present 
in the IIM model. 

In this work we perforin an analysis similar to that done in 11] by using a recently proposed parametrization for 
the dipole-proton scattering amplitude, the first developed in momentum space [l3]. 

In Section II we present a review on the main features of the QCD evolution at high energies within the dipole 
picture [lj| . Section III is devoted to describe the parametrization for the scattering amplitude in momentum space 
that will be used to describe the DIS data. The fiuctuations are properly included, the average amplitude is defined 
and the results of the fit to the last HERA data are presented. The discussion and conclusions are presented in Section 
IV. 

II. THE DIPOLE SCATTERING AMPLITUDE 
A. Dipole evolution 

Considering a dipole of transverse coordinates x and y at a given rapidity y, if one increases the rapidity from Y 
to Y + 6Y, there is a probability for a gluon, with transverse coordinate z, to be emitted by the quark (or antiquark) 
of the pair. In the large Nc limit {Nc is the number of colors), this gluon can be considered as a quark-antiquark pair 
-a new dipole- at point z. This is the dipole picture introduced by Mueller |13| . 

The probability density for the original dipole to split into the two child dipoles {x, z) and (2, y) is given by 



where a = UsNc/t^ and 



—M{x,y,z)dYdPz, (4) 

M{x^y,z)^ ^^;^^' (5) 

[x- zY{z-yY 

If the target (in this case, the proton) is dense enough, both child dipoles can interact with it and the resulting 
evolution equation for the dipole scattering amplitude is [5| 

dy {T{x, y))^aj cfzMix, y, z) [{T{x, z)) + (T{z, y)) - (T{x, y)) - (T{x, z)T{z, y))] . (6) 

This is not a closed equation for the one-dipole scattering amplitude, but the first equation of an infinite hierarchy, 
the Balitsky- JIMWLK hierarchy 0, Il3 ■ The first and second terms give the contribution of each dipole which can 
interact with the proton; the third term gives the virtual contributions to the scattering and the last suppression 
(quadratic) term comes from taking into account multiple interactions, when both child dipoles interact with the 
target. When T is small, this quadratic term can be neglected and Eq.® reduces to the (dipole version of) the linear 
BFKL equation ^. 

B. Balitsky-Kovchegov equation 

In the mean field approximation, valid when the target is sufficiently large and homogeneous, one can write 
{T {x , z)T {z , y)) w {T{x,z)) {T{z,y)) and the resulting equation is the so-called BK (Balitsky-Kovchegov) equa- 
tion |5j, J] 



dy {T{x,y)) = a J d^zM{x,y,z) [{T{x,z)) + {T{z,y)) - {T{x,y)) - {T{x,z)) {T{z,y))] . (7) 

This equation includes unitarity corrections and is free from the problem of diffusion to the infrared (nonperturbative) 
region, present in the solution of BFKL equation, since there is a scale, the saturation scale Qs(Y), an increasing 
function of rapidity, which naturally emerges from BK equation. 

If one neglects the dependence on the impact parameter, the scattering amplitude depends only on the size of the 
dipole and Eq.([7l) becomes an equation for {T{r)), where r = |r| = |cc — y\. Let us denote this mean-field amplitude 
by jVvir)- After performing the Fourier transform 



1 f d'^r f°° rlr 

iVy(fc) = ^ y — e^'^-- A/V(r) = j -Mkr)J^Y{r), 

one gets the 6-independent BK equation in momentum space 

'p^NYip)-k^NYik) k^Nrik) 



dyNy{k)^^j'4 
T^ J P 



1^2 _ p2 



which can be rewritten as 



where 



VV + fc^, 

dyNy = ax{-dL)Ny - aiV^, 

X(7) = 2V;(1) - 7^(7) - V(l - 7) 



aiV^(fc), 



(8) 



(9) 



(10) 



(11) 



is the BFKL kernel and L = log(A:^/A:o), with /cq some fixed soft scale. 

The kernel (jlip can be written with the help of a series and after a saddle point approximation and a change of 
variables fSil it has been shown that BK equation reduces to Fisher and Kolmogorov-Petrovsky-Piscounov (FKPP) 
equation [17] , which admits the so-called traveling wave solutions. In the QCD language this means that, at asymptotic 
rapidities, the scattering amplitude depends only on a single variable, the ratio k'^ /Ql{Y) instead of depending 
separately on k^ and Y . This scaling property is called geometric scaling and has been observed in the measurements 
of the proton structure function at HERA [7[ . The amplitude is a wavefront which interpolates between and 1 and 
travels towards large values of fc^ with speed A - the saturation exponent - without deformation. The position of the 
front, for which Ny = 0(1) is given by the saturation momentum Qs{Y) or, more specifically, by \og{Ql{Y) / k^) = \Y . 

The expression for the tail of the scattering amplitude reads 



Ny{k) 



k:s>Q, 



fc2 



QW) 



log 



Qi{Y) 



exp 



W {P/QliY)) 



2ax"{lc)Y 



where 



, . -X{l) -Xilc) _ ,, ^ 

A = mm a = a — ax (7c)- 



7 



7c 



(12) 



(13) 



One can observe that the last term in Eq.(|12p has an important role, since it introduces an explicit dependence on 
rapidity and hence violates geometric scaling. Then, geometric scaling is obtained for 



\og{k^lQl{Y))<^2x"{ic)aY, 
i.e., within a window Vr above the saturation scale. 



(14) 



C. Beyond the mean field approximation: the effects of fluctuations 

The BK equation is the simplest equation which describes high energy dipole evolution and scattering in perturbative 
QCD. This mean field equation has been shown to be in the universality class of FKPP equation, whose dynamics 



is called reaction-diffusion dynamics. Within this correspondence between a reaction-diffusion process and the QCD 
evolution at high energy, it has been recently realized that the Balitsky-JIMWLK hierarchy is not complete because 
they do not take into account the gluon (dipoles) number fluctuations, which are related to discreteness in the 
evolution, and thus they are completely missed by BK equation. 

As we will see in the following, the fluctuations influence dramatically the QCD evolution at high energies, and 
so the properties of the scattering amplitudes. Their inclusion results in a new hierarchy of evolution equations, the 
Pomeron loop equations i2l|. The first equation of this hierarchy is exactly the same as Eq.®, but as one goes to the 
second one, the evolution equation for the two-dipole amplitude ^T^^') = {TT), besides the linear (BFKL) term and 
the nonlinear term, responsible for unitarity corrections, there is a new term, proportional to the one-dipole amplitude 
(T), which is the fluctuation term. More generally, the equation for the fc-dipole amplitude ^T^*^^) depends on (T^^^'j, 
^j^lfc+i)^ and (T*^*^"^^), the latter being the contribution of the fluctuation. After an approximation [3| to get rid 
of the impact-parameter dependence, this can be written as a Langevin equation for the event-by-event amplitude 
which is formally the BK equation with a noise term, which lies in the same universality class of the stochastic FKPP 
equation (sFKPP). 

Each realization of the noise means a single realization of the target in the evolution and leads to an amplitude for 
a single event. Different realizations of the target lead to a dispersion of the solutions, and then in the saturation 
momentum ps = ln(Qg/A;Q) from one event to another. The saturation scale is now a random variable whose average 
value is given by 

(Q^(y))=exp[A*F]. (15) 

The dispersion in the position of the individual fronts is given by 

a^ = {pD - {ps? - DaY. (16) 

The diffusion coefficient D, as well as the average saturation exponent A*, are analytically known only in the asymp- 
totic limit a^ -^ 0, then in what follows they will be treated as free parameters. 
The probability distribution of ps is, to a good approximation, a Gaussian [l8| 



Py{Ps) - , — - exp 



{Ps-{Ps)r 



(17) 



For each single event, the evolved amplitude shows a traveling- wave pattern, which means that geometric scaling is 
preserved for each realization of the noise. However, the speed of the wave is smaller than the speed predicted by BK 
equation. This speed, or the (average) saturation exponent, has been found to be [2] 



7i'^7cX"(7c) 
ln(l/a2) 



A*^A- " '^^': - (18) 



The average amplitude is determined by {X = ln(l/r^QQ)) 



-|-oo 



(r(X,p,))-/ dp,PYips)T{X,ps). (19) 



A crucial property of the physical amplitudes is that at sufficiently high energies, unlike the individual fronts, they will 
generally not show geometric scaling. More specifically, they will show additional dependencies upon Y, through the 
front dispersion a. Then, geometric scaling is washed out and replaced by the so-called diffusive scaling |2l.l3ll4l.ll9l| 

III. DESCRIPTION OF DIS DATA 

We have seen in the Introduction that in the dipole frame the F2 proton structure function can be written in 
terms of the dipole-proton cross section adip, which can be expressed through the average dipole-proton scattering 
amplitude through Eq.Q. If one treats the proton as an homogeneous disk of radius Rp, i.e., if one neglects the 
impact parameter dependence, the amplitude depends only upon the size r of the dipole and after integrating out the 



remaining angular dependence, the dipole-proton cross section can be written in terms of the ampUtude {T(r)) in the 
following way: 



al^{r,Y)^27rRl{T{r)). 



(21) 



This expression must be inserted in Eq.([T]) and one has to parametrize the dipole scattering amplitude {T(r)) in order 
to reproduce a'^ ^ measurements. From now on we will denote this amplitude by T(r, Y). 



AGBS model and fluctuations 



The most recent parametrization for the dipole-proton scattering amplitude has been proposed by Amaral, Gay 
Ducati, Betemps and Soyez, the AGBS model IJ], which is the first parametrization in momentum space in the 
literature and is based on the knowledge of asymptotic behaviors of the solutions of BK equation. The starting 
point is that, after performing the Fourier transform ([5]), it is possible to rewrite the cross section ([1]) in terms of the 
amplitude in momentum space. The F2 structure function takes the form |12l | 



F2{x,Q^) 



47r2 



dk 
~k 



dz\^{k\z-Q')\'T{k,Y) 



(22) 



where now the photon wave function is expressed in momentum space and r(fc, Y) is the scattering amplitude in 
momentum space. The AGBS model analytically interpolates between the behaviors of the BK scattering amplitude 
in the dilute regime, which is described by Eg. p^ . and the saturation one, in which it behaves like 



f{k, Y) (k) '"'^'^ c - log 



.(Y) 



(23) 



If one defines the variable p = \n{k^ /k^), the interpolation in the AGBS model is done through the following expression 
for the scattering amplitude 



AGBS 



where 



T 



Tdii = exp 



{p,Y) = Lf (1 



-Tdii 



-lc{p- Ps) 



C'-log'm 



2ax"ilc)Y 



(24) 



(25) 



£ = ln 



1 + e'^P-P'^ 



and 



Lf = 1 + In 



with Q^(y) = fc^e 



hip~P=) _|_ g-^ip-p^) 



2 „\Y 



(26) 



(27) 



Through this parametrization, the measurements for the F2 structure function were successfully reproduced with the 
contributions of light and heavy (charm) quarks included in the fit. 

To include the fluctuations in the description of HERA data, one considers the scattering amplitude given by AGBS 
model as a single event one, and one has to evaluate the average scattering amplitude, which is obtained by performing 
the integration 



T, 



AGBS 



ipAps))) = 



+ 00 



dpsPYiPs)TP^"^ip,Ps) 



(28) 



This is the expression which must be inserted into (|22|) in order to reproduce DIS measurements of F2 structure 
function with fluctuations effects correctly included. 





xVn.o.p 


fcg (xlO-3) 


A 


R{GeY'^) 


X"(7c) 


D (xlO-2) 


rf,AGBS 


0.949 


3.79 ± 0.30 


0.213 ±0.003 


3.576 ±0.059 


4.69 ± 0.23 





^yAGBS^ 


0.949 


3.79 ± 0.30 


0.213 ±0.003 


3.576 ± 0.059 


4.69 ± 0.23 


0.0 ±1.1 



TABLE I: Parameters extracted from the fit to F2 HI and ZEUS data [2Ci.l21j in the case where mu,d.s = 50 MeV. 





xVn-o.p 


fco (xlO-^) 


A 


i?(GeV-') 


X"(7c) 


D (xlO~^) 


rf,AGBS 


0.942 


1.69 ±0.16 


0.176 ±0.004 


4.83 ±0.12 


6.43 ±0.29 





^yAGBS^ 


0.942 


1.69 ±0.16 


0.176 ±0.004 


4.83 ±0.12 


6.43 ±0.29 


0.0 ±9.6 



TABLE 11: Parameters extracted from the fit to F2 HI and ZEUS data W,^ in the case where niy, d 3 = 140 MeV. 



IV. DATA SET AND RESULTS 



In this analysis, all the last HERA data measurements of the proton structure function from HI and ZEUS CoUab- 
|20l . I2IJ are fitted, within the following kinematical range: 



orations 



X < 0.01, 



(29) 



0.045 < Q^ < 150 GeV^ 



(30) 



which corresponds to 279 data points. Both ranges include values of x low enough for the analysis to be in the high 
energy regime, and values of Q^ which allow us not to include DGLAP corrections. 

Concerning the parameters, we keep fixed a = 0.2, which enters into the amplitude through Eg. (fT3| . and jc — 
0.6275, whose value corresponds to the one obtained from the LO BFKL kernel. The other parameters in the 
amplitude, A, fcg and x", are left to be free, as well as the proton radius Rp, which fixes the normalization of the 
dipole-proton cross section with respect to the dipole-proton amplitude, and the diffusion coefficient D. Only light 
quarks are considered and the values used for their masses are iriu.d.s = 50 and 140 MeV. 

Figures [T] and [5] show the F2 structure function in bins of Q^ for small and moderate values of Q^, respectively. As 
usual, the HI data have been rescaled by a factor 1.05 which is within the normalization uncertainty. Tables HI and HIl 
show the values of the parameters obtained from the fit both with and without fluctuations. The latter corresponds 
to the value D = Q for the diffusion coefficient. 



V. CONCLUSIONS AND DISCUSSION 



In this paper, the AGBS model is used to investigate the possible effects of the gluon number fluctuations in the 
HERA data. The expression for the amplitude. Eg. ((24|) . is considered as a single-event amplitude and the average 
amplitude is evaluated through Eg. ()28p and put into the expression for the proton structure function, whose data 
was successfully reproduced. This is shown by the good x^/n.o.p. and by the curves in Figures [1] and O From the 
comparison between the results with and without fluctuations (see Tables U and |lTl) one sees that the value of the 
X^/n.o.p. does not change, and the same can be said for the parameters. Specially for the diffusion coefficient D, its 
value obtained in the case with fluctuations is very small, actually very near its mean field value D = 0. Then, in the 
framework of AGBS model, there is no evidence of fluctuations in DIS experiment at HERA energies. This indicates 
that a mean field treatment, with fixed coupling, is enough to investigate high energy QCD phenomenology, at least 
at HERA energies. 

It is interesting to compare our results with those obtained in Ref.fTli], where GBW and IIM models are used to 
parametrize the average dipole-proton scattering amplitude. Although this analysis results in a successful description 
of HERA data, it is not conclusive concerning the presence of fluctuations. In particular, the value of D extracted 
from the fit agrees with the values predicted in the literature [22, [23|, i.e., a sizable number of order C(l), which 
would indicate that fluctuations could be present at HERA. However, one should point out that in [ll|, the fit is 
performed using only data from ZEUS Gollaboration [2l|, i.e., it does not include HI data [23|, and within a more 
restricted kinematical range, with values of virtuality such that Q^ax = 50 GcV^. 

Actually, if one performs the same analysis of Ref.[ll|-only ZEUS in the same kinematical range-using the AGBS 
model, the results obtained are very similar, as it can be seen in Table IIVI with a better x^/n.o.p.. The present 
analysis is then more complete, since it includes all the last HERA data and considers a wider kinematical range. 
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FIG. 1: Predictions for the HI [20i| and ZEUS 21i| data for the proton structure function versus x for smah values of Q , given 
in GeV^. The fit was performed with quark masses mu,d,s ~ 140 MeV. 





X^/n.o.p 


fcg (xlO-3) 


A 


i?(GeV-') 


X"(7c) 


D 


jiAGBS 


0.788 


4.258 ± 0.425 


0.214 ±0.005 


3.497 ± 0.068 


4.336 ± 0.281 





/rf.AGBS\ 


0.782 


4.023 ± 0.560 


0.190 ±0.030 


3.644 ±0.214 


3.840 ± 0.214 


0.922 ±1.162 



TABLE III: Parameters extracted from the fit to F2 ZEUS data [2l| in the case where mu,d,a = 50 MeV. 





X^/n.o.p 


fcg (xlO-3) 


A 


R(GoY-^) 


X"(7c) 


D 


rf,AGBS 


0.778 


1.965 ±0.222 


0.177 ±0.006 


4.681 ±0.136 


5.946 ± 0.944 





/rf,AGBS\ 
{Ty ) 


0.768 


1.383 ±0.118 


0.120 ±0.010 


5.459 ± 0.043 


5.464 ± 0.547 


1.778 ±0.381 



TABLE IV: Parameters extracted from the fit to F2 ZEUS data [2l[ in the case where mu,d,s = 140 MeV. 

Our conclusions seem to shed some light on the investigation of the effects of fluctuations at HERA and confirm the 
robustness of the AGBS model. Of course, only in the near future, at LHC, it will be possible to see if fluctuations 
are present at much higher energies, or if they are really suppressed by the running of the coupling, as it has been 
suggested by recent developments on a one-dimensional toy model which successfully reproduces the main features of 
scattering and high energy evolution in QCD |23|.|24|. 
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FIG. 2: Predictions for the HI [20| and ZEUS [2l|] data for the proton structure function versus x for moderate values of Q^ , 
given in GeV^. The fit was performed with quark masses mu,d,s = 140 MeV. 
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